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Abstract 

In the theory of the two-dimensional Ising model, the diagonal susceptibility is equal 
to a sum involving Toeplitz determinants. In terms of a parameter k the diagonal 
susceptibility is analytic inside the unit circle, and the authors proved the conjecture 
that this function has the unit circle as a natural boundary. The symbol of the Toepltiz 
determinants was a fc-deformation of one with a single singularity on the unit circle. 
Here we extend the result, first, to deformations of a larger class of symbols with 
a single singularity on the unit circle, and then to deformations of (almost) general 
Fisher-Hartwig symbols. 


I. Introduction 


In the theory of the two-dimensional Ising model there is a quantity, depending on a 
parameter k, called the magnetic susceptibility, which is analytic inside the unit circle. It is 
an infinite sum over M, N € Z involving correlations between the spins at sites (0,0) and 
( M,N ). It was shown in [10] to be representable as a sum over n > 1 of n-dirnensional 
integrals. In [6] B. Nickel found a set of singularities of these integrals which became dense 
on the unit circle as n —>• oo. This led to the (as yet unproved) natural boundary conjecture 
that the unit circle is a natural boundary for the susceptibility. 

Subsequently [3] a simpler model was introduced, called the diagonal susceptibility, in 
which the sum of correlations was taken over the diagonal sites (A, IV). These correlations 
were equal to Toeplitz determinants, and the diagonal susceptibility was expressible in terms 
of a sum involving Toeplitz determinants. 

The Toeplitz determinant Dig(ip) is det (<Pi-j)i<i,j<N, where ipj is the jth Fourier coef¬ 
ficient of the symbol ip defined on the unit circle. The sum in question is 

OO 

Y^''PnW)-m\ 

N= 1 
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where 


¥>(0 


i-fc/e 

i-A:^ ’ 


and M, the spontaneous magnetization, is equal to (1 — A; 1 2 ) 1 / 8 . This also (as we explain 
below) is equal to a sum of n-dimensional integrals, the sum is analytic for |fc| < 1, and the 
singularities of these summands also become dense in the unit circle as n —>• oo. This led 
to a natural boundary conjecture for the diagonal susceptibility, which we proved in [9J. 

The question arises whether the occurrence of the natural boundary is a statistical 
mechanics phenomenon and/or a Toeplitz determinant phenomenon. This note shows that 
at least the latter is true. We consider here the more general class of symbols 


^) = (i-k^r + (i-k/cr~m, 


where ip is a nonzero function analytic in a neighborhood of the unit circle with winding 
number zero and geometric mean one. We assume a± Z, Rea-t < 1. The parameter k 
satisfies \k\ < 1. We define 


X( k ) = ^2 [Dn(v) ~E(tp)], 

N=1 


(1) 


where 

E(<p) = lim Lhv(<£>)0 

jV->oo 

Each summand in (JT]) is analytic in the unit disc |fc| < 1, the only singularities on the bound¬ 
ary being at k = ±1, and the series converges uniformly on compact subsets. Therefore 
x(k) is analytic in the unit disc. 

Theorem. The unit circle \k\ = 1 is a natural boundary for x(k). 

The result in [9j was established, and here will be established, by showing that the 
singularities of the nth summand of the series are not cancelled by the infinitely many 
remaining terms of the seriesH We shall see that a certain derivative of the nth term is 
unbounded as k 2 tends to an nth root of unity while the same derivative of the sum of the 
later terms is bounded, and if it is a primitive nth root the same derivative of each earlier 
term is also bounded. 

To put what we have done into some perspective, we start with a symbol 

(1-0“ + (1-1/0°-• 

1 When ?/>(£) = 1 it equals (1 — k 2 )~ a + a ~. In general it equals this times a function that extends 
analytically beyond the unit disc. 

2 A nice example [7] where such a cancellation does occur is 

z _ z z 2 z 4 z 2 " 
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Then we introduce its ^-deformation, times a “nice” function and consider their 

Toeplitz determinants as functions of the paremeter k inside the unit circle. Is it important 
that we begin with a symbol with only one singularity on the boundary? It is not. We may 
begin instead with a general Fisher-Hartwig symbol [2j 

P Q 

n a-",*)”* ii(i 

P= 1 9=1 

where \u p \, \v q \ = 1 and P, Q > 0. With some conditions imposed on the a+ and the a~, 
we show that the conclusion of the theorem holds for the deformations of these symbols. 

Here is an outline of the paper. In the next section we derive the expansion for y(/c) 
as a series of multiple integrals. In the following section the theorem is proved, and in 
the section after that we show how to extend the result to (almost) general Fisher-Hartwig 
symbols. In two appendices we give the proof of a proposition used in Section II and proved 
in [5], and discuss a minimum question that arises in Section IV. 

II. Preliminaries 


We invoke the formula of Geronimo-Case [5j and Borodin-Okounkov [2j to write the 
Toeplitz determinant in terms of the Fredholm determinant of a product of Hankel opera¬ 
tors. The Hankel operator is the operator on £ 2 (Z + ) with kernel ((pi + j + N+i)ij>o- 

We have a factorization y>(£) = <p + (£) </?_ (£), where <p+ extends analytically inside the 
unit circle and </?_ outside, and <^ + (0) = </?_(oo) = 1. More explicitly, 


<p + (x) = (l-kO a+ ^ + (Z) and <p_(0 = (l-k/0 a ~i>- CO- 


If V>(£) is analytic and nonzero for s < |£| < s 1 then ip+(£) resp. 4>~(C) is analytic and 
nonzero for |£| < s^ 1 resp. |£| > s. 

The formula of G-C/B-0 is 


Dn{<p) = E(<p) det (j - ^A r (^ 1 )) 


where for a function / we define /(£) = /(£ x ). Thus, if we write 

A(f) = = C 1 - k£)- a + (1 - k/£) a ~ 


<P+(Z) 

Kn = Hn(A) Hn( A -1 ) 


V’+W’ 


( 2 ) 


then x(k) equals E(ip) times 


OO 
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In |9J the following was proved. We give the proof in Appendix A. 

Proposition. Let Hjy(du) and H^{dv) be two Hankel matrices acting on £ 2 (Z + ) with i,j 
entries 

J x N+i+ > du(x), j y N+i +> dv(y), (3) 

respectively, where u and v are measures supported inside the unit circle. Set Kn = 
H^{du) H]y(dv). Then 


^[det(/-A»-l] 

JV=1 


\ (~l) n [ ... [ Ilj x iVi 

( nl ) 2 J J i-FI i x iVi 


det 


1 


1 - Xiyj 


2 

J\du(xi) dv(yi), 
i 


where indices in the integrand run from 1 to n. 


We apply this to the operator Kn = Hat (A) Hjy(A *) given by ©. The matrix for 
Hn (A) has i.j entry 

T / a or—* 

where the integration is over the unit circle. The integration may be taken over a circle 
with radius in (1, as long as \k\ > s. (Recall that r ip± (£) are analytic and nonzero for 

s < |£| < s^ 1 .) We assume this henceforth. 

Setting ( = l/i we see that the entries of Hn( A) are given as in Q with 


du(x ) =- A(x x ) dx, 

2m 


and integration is over a circle C with radius in (|/c|, 1). Similarly, Hn(A 1 ) = Hn(v) where 
in © 

dv{y) = -—: A(y) _1 dy, 

2m 

with integration over the same circle C. 

Hence the Proposition gives 

OO 

S(k) = ^Sn(k), (4) 

71=1 

where 


S n (k) 


(-If 1 

(n!) 2 (2iri) 2n 


(det ( _ — 

1 - rii x iVi V V 1 - x iVj 


2 TT ^ Xi ^ 

11 Mvi) 


Q dx -, dijj . 
i 


with all integrations over C. 
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We deform each C to the circle with radius |/c| (after which there are integrable sin¬ 
gularities on the contours). Then we make the substitutions Xi —>• kxi, yi —>• kyi, and 
obtain 


S n (k) 


(— l) n K 2n 
(n!) 2 (27ri) 2n 


, 11 (det (-—1- 

V V 1 ~ KX iUj 


2 TT Mk lx i *) 

J-.J- A (ky^ 


where integrations are on the unit circle. We record that 

A (k~ 1 x~ 1 ) (1 — Kx) a ~ (1 — x~ l )~ a + p(fc -1 a; -1 ) 

A (ky) (1 - Ky)~ a + (l-y _1 ) Q - p(ky) 


Y[dxi dy,, 

i 

(5) 

( 6 ) 


where we have set 


k = k 2 , p(x) 


V'-Qk) 

^+(x)' 


The complex planes are cut from k 1 to oo for the first quotient in <[6|) and from 0 to 1 for 
the second quotient. 


Using the fact that the determinant in the integrand is a Cauchy determinant we obtain 
the alternative expression 


9 iiA = (~ 1 ) n Kn{n+1) f ... f Eh xm A(x) 2 A (y) 2 

(n!) 2 (2vr i) 2n J J 1- K n T\i x iyi Yliji 1 - KXiyj) 2 

where A(x) and A(y) are Vandermonde determinants. 


TT Mk 1 x i 
J-.J- A (ky^ 


dxidy 


i i 


i 


(7) 


For any 5 < s we can deform each contour of integration to one that goes back and forth 
along the segment [1 — S, 1] and then around the circle with center zero and radius 1 — 
This is the contour we use from now on. 


III. Proof of the Theorem 


There will be three lemmas. In these, e ^ 1 will be an nth root of unity and we consider 
the behavior of S{k ) asfc->e radially. Because the argument that follows involves only the 
local behavior of S(k), we may consider k as the underlying variable and in ([6|) replace k 
by the appropriate yfn. We define 

[i = K~ n — 1, f3 = a + + ot-, b = Re/3, 


so that p > 0 and p —>• 0 as k —>• e. 


Lemma 1. We have@ 



2n 2 — [bn] 

Sn(k) 


[bn]-Pn-l 


3 To expand on this, it goes from 1 — <5 to 1 just below the interval [1 — 8, 1], then from 1 to 1 — <5 just 
above the interval [1,1 — <5], then counterclockwise around the circle with radius 1 — J back to 1 — J. 

4 We use the usual notation [bn] for the greatest integer in bn. The symbol « here indicates that the ratio 
tends to a nonzero constant as y, —> 0. 
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Proof. We set 

l = 2 n 2 — [bn] 

and first consider 


Ui x m 


A(x) 2 A (yf 


(1 - K n n f XiyiY +l n M -(l - «®i3/j) 2 


n 


A(k 


X X: 


A {kyi 


j I dxi dyi 


( 8 ) 


where all indices run from 1 to n. This will be the main contribution to d^S n (k)/diY. 

For the i,j factor in the denominator in the second factor, if Xi or y.j is on the circular 
part of the contour then \x % y,j <1 — 5 and the factor is bounded away from zero; otherwise 
Xiyj is real and positive and this factor is bounded away from zero asK->6 since e / 1. So 
we consider the rest of the integrand. 

ifn < \xiy%\ <1 — 5 then the rest of the integrand is bounded except for the last quotient, 
and the integral of that is 0(1) since Rea± < 1. 

When f](- \xiyi\ >1 — 5 then each |xj|, \yi\ > 1 — 5, so each Xi, yi is integrated below and 
above the interval [1 — 5,1]. If all the integrals are taken over the interval itself we must 
multiply the result by the nonzero constant (4 sin7ra + sinTra-)" - . The factors 1 — nxiyj in 
the second denominator equal 1 — /t(l + 0(5)) = (1 — k) (1 + 0(5)) since k is bounded away 
from 1. From this we see that if we factor out n( i+1 ') 71 from the first denominator, (1 — n) n 
from the second denominator, and (1 — Ac)^ n ( / o(fc~ 1 )/ / o(/c)) ?1 from the last factor (all of these 
having nonzero limits as n —>• e), the integrand becomes 

no. - *)-(i - + ow). 


We make the substitutions Xi = 1 — £j, yi = 1 — ry and set r = JT(£j + yi). Then since 
n f (l — ^i)(l — rji ) = 1 — r + 0(r 2 ) this becomes 


A(Q 2 A(, ? ) 2 

{n + r + 0(?’ 2 ))^ +1 


(i + o(5)). 


The integration domain becomes r < 5 + 0(5 2 ). Consider first the integral without the 
0(5) term. By homogeneity of the Vandermondes and the product, the integral equals a 
nonzero constant times 


L 


8+0{5 2 ) 


r 2n 2 —(3n—l 


(y + r + 0(r 2 )y +1 


dr. 


(9) 


Making the substitution r —>• yr results in 


5 This is the integral of A(£) 2 A(r;) 2 ]”[ ^ a+ rj i a over r = 1. It can be evaluated using a Selberg integral, 
with the result 


1 

r(2n 2 - /3n) 


Il r (i + 2) 2 r(j 

3=0 


a+ + 1) r (j - a- + 1). 
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( 10 ) 




2 n 2 —/3n—£— 


'L 


(5+0(«5 2 ))/ M 


r 2n 2 —f3n— 1 


(1 + r + 0(p 2 r 2 )Y +1 


dr 




‘1 


(5+0(5 2 ))//, 


r ,2n 2 —/3n—1 


(1 + r + 0(^2 r 2))2n 2 -[6n]+l 


dr, 


where we have put in our value of t. 
integral 



The integral has the p —» 0 limit the convergent 
2n 2 —[im]+l 

| r [&n]-/2ra-2 ^ 


and © is asymptotically this times pY™^ 1 . 

For the integral with the 0(6) we take the absolute values inside the integrals and 
find that it is 0(6) times what we had before, except that the (3 in the exponents are 
replaced by 6, and in footnote [5] the exponents a± are replace by their real parts. Since <5 
is arbitrarily small, it follows that the intergral of ([9]) is asymptotically a nonzero constant 
times 

To compute the derivative of order 2 n 2 — [bn] of the integral in (J71) one integral we get 
is what we just computed. The other integrals are similar but in each the t in the first 
denominator is at most 2 n 2 — [bn] — 1 , while we get extra factors obtained by differentiating 
the rest of the integrand for S n (k). These factors are of the form (1 — KXiyi ) , (1 — KXi ) -1 , 
(1 — Kyi ) -1 , or derivatives of p(k~ 1 x^ 1 ) or of p(kyi)~ 1 . These are all bounded. Because 
£ < 2 n 2 — [bn] — 1 the integral (|9|) is 0(p~ l+1 ) for some 7 > 0. The lemma follows. □ 


Lemma 2. If e m / 1 then 


7 \ 2 n 2 — [bn\ 

e) «■<*)= °(D- 


Proof. If e m Y 1 terms, aside from those coming from the last factors, obtained by 
differentiating the integrand in ([7]) with n replaced by m are bounded as k —» e. Differenti¬ 
ating the last factor in the integrand any number of times results in an intregrable function. 
□ 

Lemma 3. We have 

( 7 \ 2 n 2 — [bn] 

*) «■<*> = °w- 

Proof. We shall show that for k sufficiently close to e all integrals we get by differentiating 
the integral for S m (k) are at most A m m m , where A is some constant]! Because of the 
l/(m !) 2 appearing in front of the integrals this will show that the sum is bounded. 

6 The value of A will change with each of its appearances. It may depend on n and 5, which are fixed, 
but not on m. 
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As before, we first use (J7|) with n replaced by m, and consider the integral we get when 
the first factor in the integrand is differentiated 2ra 2 — [bn] times. All indices in the integrands 
now run from 1 to m. 


First, 

|1 - n m \\_ x iVi\ > 1 - 

i i 


Next we use that either \xi\ = 1 — 5 or x% € [0,1], and n € [0, e], to see that |1 — KXi\ > 
min(5, d), where d = dist(1, [0, e]). We may assume 5 < d. Then \l — KXi\ > 5 , and similiarly, 
|1 — nyi\ > 5. It follows that the integrand in dTJ) after differentiating the first factor has 
absolute value at most A m times 


1 


A(x ) 2 A {yf 


(i-Ili\xiyi\) 2n2 - [bn]+1 11./ Ii 

where a± = Re a±. 


nn ~ Xi \ a+ 1 1 ~ ~ 


(ii) 


If n,- \xiDi\ <1 — 5 then the hrst factor is at most When J~[ t \xiVi\ >1 — 5 

we set, as before, Xi = 1 — £*, yt = 1 — rji with rji € [0, 5]. Since we are to integrate back 
and forth over these intervals we must multiply the estimate below by the irrelevant factor 

22m 


We have ]/[,(! — &)(! — Vi) < (1 — £i)(l — Vi) for each i, and so averaging gives 
lid - &)(1 - Vi) < ^ 1 ~ &)(1 - Vi), 

i i 

and therefore 

1 - IF 1 - € 0(1 - w) > A - (1 - € 0(1 - n)) 

i i 

= 7^ + Vi- ZiVi) > 2^ 5^(6 + ^)/ 2 ( 12 ) 

i i 

if 5 < 1/2, since each < 5 /, r/j < 5. From this we see that in the region where + Vi) > 5 

the first factor in (HID is at most (4 m/5) 2n2 [M+ 1 . 

So in either of these two regions the first factor is at most A m . We then use (HID with 
the second factor replaced by the absolute value of 




Each denominator has absolute value at least 5, so by the Hadamard inequality the square 
of the determinant has absolute value at most 5~ 2m m m . Therefore the integral over this 
region has absolute value at most 


A m m n 


J ■■■JH\l-x i \- a +\l-y i \- a ~ 


0 dxid yi . 

I 





The integral here is A m , and so we have shown that the integral in the described region is 
at most 

It remains to bound the integral over the region where Xi = 1 — Vi = 1 — rji with 
£ i. r)i € [0,5], and r = £T(£j + Vi) < <5- Using (fT2J) again, we see that the integrand has 
absolute value at most A m times 

d -m 2 A(Q 2 A(r?) 2 "1 T £-a + -o_ 

(E,:te + ^)) 2n2 - [H+lii ^ Vi ’ 


(Recall that d = dist(1, [0, e]), snd KXiUj € [0, e]. The factor (4m 2 ) 2 " 2 coming from 

using m were absorbed into A m .) Integrating this with respect to r over r < 5, using 
homogeneity, gives 

A (0 2 A iv) 2 II^~ a+ ^~ a " d ^’V) 

1 i 

(where d(£, rj) denotes the ( 2 n — l)-dimensional measure on r = 1 ) times 



d ~ m2 



2m 2 —In 2 + [bn]—bn— 1 


dr. 


The first integral is given in footnote [5] with n replaced by m and a± replaced by a±, and 
is exponentially small in m. The last integral is 0(S 2m ") since m > n and n is fixed. Since 
5 2 < d, the product is exponentially small in m. 

So we have obtained a bound for one term we get when we differentiate 2 n 2 — [ 7 n] times 
the integrand for S m (k). The number of factors in the integrand involving k is 0(m 2 ) so if 
we differentiate 2 n 2 — 1 times we get a sum of 0(m 4n ) terms. In each of the other terms 
the denominator in the first factor has a power even less than 2 n 2 — [ 771 ] and at most 2 n 2 
extra factors appear which are of the form (1 — nxiyi)^ 1 , (1 — KXj) -1 , or (1 — Also, 

^(A; -1 ^^ 1 ) or p(kyi )~ 1 may be replaced by some of its derivatives. Each has absolute value 
at most 5 _1 , so their product is 0(5 -4n2 ). It follows that we have the bound A m m m for 
the sum of these integrals. Lemma 4 is established. □ 


Proof of the Theorem. Let e be a primitive nth root of unity. Then e m 7 ^ 1 when m < n 
so Lemma 2 applies for these m. Combining this with Lemmas 1 and 3 we obtain 


^ \ 2n 2 — [bn] 
dK 


s(k ) « 1 


as k —>• e. This is unbounded, so S(k) cannot be analytically continued beyond any such e, 
and these are dense in the unit circle. 

Thus the unit circle is a natural boundary for S(k), and this implies that the same is 
true of %( k). □ 


9 



IV. Fisher-Hartwig symbols 


In this section we show how to extend the proof of the theorem to deformations of 
Fisher-Hartwig symbols. 

We start with a Fisher-Hartwig symbol] 

P Q 

n na-vo 0 ', 

P= 1 9=1 

where \u p \, \v q \ = 1 and P, Q > 0, and then its /c-deformation 

P Q 

no =n - ku p® at na- wo a * • 

P =1 9=1 

We assume that Rea+, Rea~ < 1 and a+, a“ 0 Z. (Plus a simplifying assumption that 
comes later.) 

The singularities of Dn(lp) on the unit circle are at the {u p v q )~ 1 / 2 , and 

E (<p)=H(i - k 2 Up v q )- a >; . 
p,q 


We have now 


HO = lid - ku P 0- a * (1 - kv q /O a \ 

p,q 


A (k l x x ) 
A (ky) 


y-r (1 — U p /x ) a P (1 — KV q x) a 1 

p,q (1 - Ku p y)-H (1 - v q /y) a « 


Again we begin by considering the integral 

f f n, x iVi A(x) 2 A(y) 2 


(1 - K n [I, XiyiY +l n M (! - 


n 


-i -i> 


A (k 1 x i 
A (kyi 


11 <ixi (bu¬ 


rn 


Our integrations are for the Xi around the cuts [1 — <5,1] u p and for the y* around the 
cuts [1 — 6, 1] v q and then both around the circle with radius 1 — 6. (In case we do want to 
generalize with a factor ip(0 as before.) If we replace integrals around the cuts by integrals 
on the cuts, then for a cut [1 — 6, 1] u p we must multiply by 2 sin 7rand for a cut [1 — 6, 1] v q 
we multiply by 2 sin7ra“. (These are both nonzero.) We assume that this has been done. 

We now let n —y a radially, where e is an nth root of n( u Pi %) > but not equal to 

any {u p v q )~ l . We also choose it so that it is not an mth root of any product of the form 
D[ (u Pi VqH 1 with m < n . These e become dense on the unit circle as n —» oo. The last 

'We could easily add a factor ip(£) to give the general Fisher-Hartwig symbol. 
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condition assures that the integrals with m < n are bounded, which will give the analogue 
of Lemma 2. We now consider the analogue of Lemma 1. 

The integral over \xiyi\ < 1 — <5 is bounded, as before. In the region where n \xiyt\ > 
1 — 5 each Xi and y^ is integrated on the union of its associated cuts. This is the sum of 
integrals in each of which each Xi is integrated over one of the cuts and each yi is integrated 
over of the cuts. Suppose that Xi is integrated over [1 — 5, 1] u Pi and y* is integrated over 
[1 — 6, 1] v qi . (We consider this one possibility at first. Then we will have to sum over all 
possibilities.) 

If we factor out ]~[ u ih from the first numerator, [7(1 — KU Pi v qi ) 2 from the second 
denominator, and na — KUp i Vq i )°‘ Pi+api from the last product the integrand becomes 
1 + 0(5) times 

(1 n (1 - u ° Jx ' r< (1 - •*'*>-*■ (14) 


We make the substitutions X { = (1 — £j) u Pi , y* = (1 — ry) v qi , and define 

Ip = {i‘-Pi=p}, Iq = {i ■ Qi = q}- 

Then 

n (i-u Pi / Xi )-<(i-v qi / yi )-^= n er 4 - n +om 

i P,i£lp <?,*€+ 


As for the Vandermondes, we have 

A(x) = ± JJ A ( Xi : i & I p ) ■ U JJ (xj - xy), (15) 

p p^p' jei P ,j'ei p r,j<j' 

and similarly for A(y). If we define n p = \I p \, n q = \I q \, then the last double product is to 
within a factor 1 + 0(5) equal to 


± n k-v) w > 

p<p' 

while the first product is to within a factor 1 + 0(5) equal to 

J] np p(np+1)/2 J1 A (& : i € I p ). 

P p 


Thus, if we factor out (re n [7 u Pi v qi Y +1 from the denominator in cm, and set y = 
K~ n Uiu Pt V qi) 1 ~ 1, then Cm nmy get replaced by a constant times 1 + 0(5) times 


1 

(M + Eite + %)) m 


A (b-.ielp) 2 ^ 0 * 

P,i£lp 


J] A(yj : i € I q ) 2 rji °^. 
q,i£lq 


(16) 
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If we use homogeneity the integral of (1161) becomes a nonzero constant ti 


times 


f 


,- 1 + J2p n p( n p- a t) + J2 q n q( n q- a q ) 


/o (/x + r)^ 1 

This is largest when the power of r is smallest. So we minimize 

Y n p {n p - dp), (dp = Re a+) 


( 17 ) 


over all { n p } with n p > 0, Yl p n p = n - The solution are not necessarily unique. But in any 


case 


M n ■= min Y n P (n p - a+) = — + 0(n) 


(18) 


and MT j > M+ for large enough n@ 
Similarly, with a~ = Re a~ and 


M n := min Y n q( n g ~ a q) = -Q+ °( n )- 


(19) 


Then we choose 

£ = J2 n P + J2 n q~ [52 n P a p + X/ ri( i a Q 

p q p q 


with the minimal n p and n q . The integral (1171) is equal to 

+ n P a t + E 9 n q a q}-( Ep n P a t +Eq n q <*q ) 


( 20 ) 


times 


/ 




P 1 + Ep n p( n p a j()+Eq n 5( n 9 OLq ^ dr 


I o (l + 0 m 

The exponent of (i has real part in (-2,-1] and the integral has a nonzero limit (a Beta 
function) as fi —» 0. 


Once we take care of the integrals with the 0(5) as in the proof of Lemma 1 we deduce 
that this is the asymptotic result for the integral when we choose this set of cuts. 

We now assume the minimal solutions are unique ^ I 


8 Also computable using a Selberg integral, it is 


_1_ 

r (E p n p( n p ~ a t) + E„ n q {n q 


a p )) 


rip — 1 riq — 1 

n n r ^'+ 2 ) r ^' - < 4 + 1 ) n n r o-+ 2 ) r u - + *)■ 

p j=0 q 3=0 


9 See Appdendix B. 

10 We shall see in Appendix B that for large n uniqueness is a condition on the <Zp and a q that depends 
only on the residue classes of n modulo P and Q. It suffices for our purposes that we have uniqueness for 
some sequence n —» oo. 


12 







Then for the other choices of cuts the integral (fT71) is 0(p~ 1+1 ) for some 7 > 0, and so 
the integral over the chosen set of cuts dominates. We still have to allocate the x\ and y % to 
the various cuts, once the numbers of each have been chosen. The number of ways of doing 
this is n!/J]n p ! for the Xi and n\/\\n q \ for the y{. (The total number of ways is at most 
/’" (,)".) 

This takes care of the integral (ITSD . the main contributions to (d/duYS n (n). We com¬ 
plete the proof of the analogue of Lemma 1 as we did at the end of the proof of that lemma. 
Thus, with t given by (f20l) . 

\t 

) S n (k) fj J ~ 1+ ^P n P a P +T, q nq a q ]-(T, p n P a++J2 q n q a q ) _ 



For the analogue of Lemma 3 we first consider the integral (1131) with n replaced by m > n, 
and £ given by (120j) . As before it remains to bound the integrals over the regions where 
each Xi — (1 £i)u Pi and each y i = ( 1 - r]i)v qi , with &, ip € [0, <5], and r = + Vi) < <*>■ 

Replacing the first denominator in (fl3l) by (YY&+Vi)Y +1 introduces a factor (4 m 2 /dY +1 
as before, a factor that can be ignored. The reciprocal of the second denominator is at most 
d~ m where d = min Pi9 dist([ 0 , e], (u p v q )~ 1 )- The product of the terms involving k in the 
last product is d~°^ m \ and so may also be ignored. The square of the product over p < p' 
in (USD, times the square of the analogous product over q < q' , is at most 2( p 2 +( 3 2 ) m2 . There 
remains an integrand whose absolute value is bounded by 


1 

+ Vi)Y +1 


n A 

p,i&I P 


A(r h :iel q ) 2 r h 0 -. 

q, iGlq 


The integral of the products over r = 1 (given exactly in footnote [HD is trivially at most 

2 

its maximum (at most A m 4 m ) times the (2m — 1 (-dimensional measure of r = 1, which is 

2 

1/T(2m). We use the crude bound 4 m . This is to multiply 

r -^-2+E p m p (m p -a+)+J 2 q m q (m q -a q ) ^ 

Now 

^m p (m p - a p ) + X/ m q (rri q - a~) 
p q 

is at least M p + M~ , and it follows from (fl8l) and (fT9l) , and the strict monotonicity of the 
sequences {M±}, that for large enough n and some R this greater than i + 1 + m 2 /R for 
all m > n. Then the integral is at most 5 m ^ R . 

This integral is one of at most P m Q m integrals, and this factor also can be ignored. 
The factors we had before that could not be ignored combine to (4 2 p2+ Q 2 /d)™ 2 . It follows 
that if we choose 5 < (4 2 p +< ^ /d)~ R the integral over r < 5 of (1131) with m replacing n is 
exponentially small. 

This takes care of the integral (ED with m replacing n, the main contributions to 
{d/dnYS m (K). We complete the proof of the analogue of Lemma 3 as we did at the end of 
the proof of that lemma. □ 
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Appendix A. Proof of the Proposition 


The Fredholm expansion is 


00 /_ 1 \n 

det (/ — Kjsr ) = 1 + EE^ E det(K N (pi,pj)). 


n\ 

n= 1 pi,...,p n >0 


Therefore its suffices to show that 


E E det(K N (p il p j )) 

N=l pi,...,p„>0 


1 

n\ 


n ‘ xm (&*,< 1 


1 - n x iVi V V 1 - 


dn(xi) • • • dn(x n ) dv(yi) ■ ■ ■ dv{y n ). 


We have 


X N+Pi yN+Pj 

Kn{PuPj)= I I - _ — du(x) dv(y). 


It follows by a general identity [Tj (eqn. (1.3) in [8j) that 

det (K N (pi,pj)) = — f ••• f det(xf +Pj ) det (yf +Pj ) \ du{xi) dv{yi) 

Til J J . i XiPi 

1 f 


n\ 


n^f det(x * j) det ^ j) n n du(xi) dv(yi 

' J- XiPi 


Summing over N gives 


^2 det (K N (p i ,p j )) 


N =1 


h.J-J iwrfrlk det<a: ‘ i) d * ( ^> 

(Interchanging the sum with the integral is justified since the supports of u and v are inside 
the unit circle.) 

Now we sum over pi,... ,p n > 0. Using the general identity again (but in the other 
direction) gives 


det(rc^) det(yf d ) = n\ det J x p p p - J = n\ det ( 

pi,...,p„>0 \p> 0 / ^ 


1 


1 - XiPj 


We almost obtained the desired result. It remain to show that 

1 \ XT 1 


det 


n 


1 - XiPj J 11 1 - XiPi ’ 


( 21 ) 
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which we obtain in the integrand, may be replaced by 


< 22 > 

This follows by symmetrization over the x,;. (The rest of the integrand is symmetric.) For 
a permutation tv, replacing the x t by multiplies the determinant in (|2T1) by sgn7r, so 
to symmetrize we replace the other factor by 


1 

n! 


£sgn- 

7r 


1 


1 •E'ir(i)yi 


— det 
n! 



Thus, symmetrizing (1211) gives (1221) . 


□ 


Appendix B. The minimum question 


Changing notation, we consider 

k 


= min 


{ ^2 rii (m - ai ) : n % E Z + , y n, = n j. 


2 = 1 


2 = 1 


and ask when this is uniquely attained. Set 

k 


s = k 1 y a*, ai = (ai — s)/2 , m = rii - n/k , 


and define 


2=1 


N fc = 


K 

{( Xi )£R k :Y,Xi = 0 }. 


2=1 


Then a = (cb) € and n = (hj) E N fc . If n = zz (rnod/c) the other conditions on the h, 
become 

hj > —n/k, rii € Z — n/k. 

(Think of z/ as fixed and n as variable.) A little algebra gives 


y rii (rii - ai) = y(n* - a*) 2 + fc (ra/fc - s/2) 2 - y a 2 /4. 


2=1 


2=1 


2=1 


Minimizing the sum on the left is the same as minimizing the first sum on the right, 
with the stated conditions on the fq. Several things follow from this. First, since the 
minimum of the first sum on the right is clearly 0(1), the condition fq > —n/k may be 
dropped when n is sufficiently large; second, M n = n 2 /k — sn + 0(1); third (from this), 
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M n+ \ — M n = 2 n/k + 0(1) > 0 for sufficiently large n; and fourth, for uniqueness we may 
replace our minimim problem by 

k 

min | ^^(rh — (h) 2 : n € N fc , hi € Z — z//&j. 

i=i 

This minimum is uniquely attained if and only if there is a unique point closest to a in the 
set of lattice points (Z — v/k) k in N fc . This condition depends only on the residue class of 
n modulo k. 

When k = 2 the subspace N 2 is the line x\ + x^ = 0 in R 2 . When n is even the 
lattice consists of the points on the line with coordinates in Z and a is equidistant from 
two adjacent ones when ai — 02 € 4Z + 2; when n is odd the lattice consists of the points 
of the line with coordinates in Z + 1/2 and a is equidistant from two adjacent ones when 
ai — a 2 € 4Z. Non-uniqueness occurs in these cases. 
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